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Summary. — I give an overview of the present knowledge about nonperturbative 
functions parametrizing the fragmentation into one or two hadrons of (un) polarized 
light quarks in vacuum, including information on their transverse momentum de- 
pendence. 

PACS 13. 66. Be - Hadron production in e'^e~ interactions. 

PACS 13 . 87 . Fh - Fragmentation into hadrons. 

PACS 13 . 88 . +e - Polarization in interactions and scattering. 



1. — Introduction 

The fragmentation process describes the transition from a highly virtual parton i at a 
scale to one hadron h carrying a fraction z of its energy. The information is encoded 
in the fragmentation function Di^{z,Q'^), which is a nonperturbative object since the 
hadronization itself is a soft, nonperturbative, process. As such, _D^j(z, Q^) cannot be 
deduced from first principles but must be extracted from experiments. 

A large amount of data has been collected in the last 30 years by measuring hadron 
spectra in e+e^ annihilations. Based on these data, several parametrizations of (z, Q^) 
have been released. More recently, new measurements in Semi-Inclusive Deep-Inelastic 
Scattering (SIDIS) and in hadronic {p-p and p-p) collisions were included in various fits. 
A very brief overview is given in Sec. 2 and the most updated parametrizations of D^^ 
are compared in Sec. 2T (for brevity, only for light partons in vacuum; for a short review 
on medium modifications and heavier flavors, see Ref. [1]). 

The dependence of ^ upon the transverse momentum fcy of the fragmenting parton 
is basically unknown. Most of the phenomelogical studies are based on a simple flavor- 
and ^;-independent Gaussian ansatz. But in several experimental results for hadron 
multiplicities the evidence emerges about transverse-momentum distributions depending 
on both energy and flavor of the detected hadron. This topic is directly addressed in 
Sec. 2'2, and also in Sec. 3 where the main three types of models of fragmentation 
functions arc sketched. 

The above considerations apply also to polarized fragmentation functions, actually 
to the only one that has been parametrized so far: the Collins function. In fact, only 
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its first kr moment could be extracted leaving the kr dependence fully unconstrained 
(see Sec. 2 '3). The Collins effect in spin asymmetries in SIDIS is one crucial tool to 
address the socallcd transvcrsity parton distribution [2], a poorly known cornerstone in 
the knowledge of the (spin) partonic structure of the nucleon. That is why on one side 
models of the Collins function were developed and studied in detail (see Sec. 3'1 and 3'3), 
and on the other side alternatives were considered in the " hunting for transvcrsity" . 

The most promising alternative is based on a spin asymmetry in SIDIS with two 
hadrons detected in the final jet. The corresponding Di-hadron Fragmentation Functions 
(DiFF) arc encoded in functions like D'1\'^^{zi, Z2, M^, Q"^), that must depend also on the 
invariant mass of the hadron pair, ; the latter represents a second natural scale in the 
fragmentation, with <C [3]. Further details and some first results are presented 
in Sec. 4. 

2. — Single-hadron Pragmentation Functions 

In order to extract information on D'l^{z,Q'^) from data, the most suitable process 
is by far the electron-positron annihilation. Measuring the socalled scaled-energy distri- 
bution (l/ctot)^- gives direct access at leading order (LO) in as to the fragmentation 
function summed over all active flavors [4]. Well established factorization theorems [5] 
allow to explore higher orders in terms of perturbatively calculable coefficient functions, 
that are known up to NNLO in the MS scheme [6, 7]. A large amount of data has 
been collected in the last 30 years in a wide energy range, 12 < Q < 200 GeV and 
0.005 < z < 0.8, and for various hadron species: tt^, 7^=*=, _ftr".p,p, A, A (see Ref. [8] 
and references therein, for a short review). Most experiments wcxv, able to disentangle 
the contribution of light quarks (u, d, s) from c and b jets. In particular, the OPAL 
collaboration released also an analysis with full flavor separation [9]. 

However, since at LO the e^er annihilation leads to the back-to-back production of 
a quark and an antiquark jet, data only allow for the extraction of the flavor-inclusive 
fragmentation function D^^ -\- £>f^(= D^ ,^ + Di ^). Moreover, the gluon fragmentation 
function D'l ^ can only be extracted from 3-jet events that by construction appear at 
NLO; hence, it is weakly constrained. 

Fortunately, these drawbacks can be compensated by Semi-Inclusive Deep-Inelastic 
Scattering (SIDIS) data as well as by data on hadronic collisions. In the valence region 
(xb > 0.1), quarks are produced more abundantly than antiquarks and individual D^^^ 
or Di g can be independently extracted. In hadroproduction with p-p and p-p collisions, 
the Di g can be directly addressed, particularly for xb <^ ^ or, equivalently, for hadron 
transverse momenta P± small with respect to the center-of-mass (cm) energy y/s available 
in the collision. Moreover, since xb = 0{P±/z) the much larger parton densities of 
the projectile/target at small xb allow to probe the fragmentation functions at large z 
(> 0.7), complementing the information extracted in e+e~ annihilations. 

SIDIS data have been collected in the last 15 years for both unidentified (h^) and 
identified charged hadrons {Tr^,K^, and also A, A) mostly in e~p collisions at HERA 
(HI [10, 11, 12], HERMES [13], and ZEUS [14, 15, 16] collaborations), also at CERN 
with muonic (anti)neutrino beams (NOMAD [17]). The explored kinematical range, 
1 < Q < 100 GeV and 0.1 < 2 < 1, signiflcantly enlarges the phase space available 
to D'i ^{z,Q'^), since the hard scale is not constrained at the cm energy as in e+e" 
annihilations, Q = \/s/2. By analyzing scaled-energy distributions in the Breit frame 
for events in the kinematical current region at z = Ph/iQ/'i), it was possible to compare 
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Table I. - Main features of global fit analyses DSS [27, 28, 29], HKNS [26], and AKK08 ]30]: 
data sample, kinematic range covered, technique for error analysis. 



DSS 


HKNS 


AKK08 




c+e", SIDIS, pp 
0.05 < 2 , 1 < < 10^ GeV^ 


0.01 < z , 1 < < 10* GcV^ 


e+e , pp, pp 

0.05 < 2 , 2 < Q2 < 4 X lO"* GeV' 


Lagrange multipliers 


Hessian errors 


in progress 



the results with the corresponding e+e~ ones and to successfully test the universality of 
fragmentation functions [12]. 

Hadron spectra in hadronic collisions appeared more recently, thanks to high-precision 
p-p measurements at RHIC (BRAHMS [18], PHENIX [19, 20], and STAR [21, 22, 23] 
collaborations) and p-p ones by CDF [24, 25] at the Tevatron. At RHIC, P± distributions 
of tt's, JT's, and protons, were measured up to 10 GeV at mid-to-large rapidities (for n^, 
also up to 20 GeV by PHENIX [20]), as well as puzzling data for A, A, by the STAR 
collaboration [23]. 

2'1. Unpolarized fragmentation. - The year 2007 represents a sort of tmning point 
for the phenomenological work about extraction of D^]^ ^{z,Q^) from experiments. All 
parametrizations released before this date are based on e+e~ data only, they suffer from 
large uncertainties at large z and Q^. and fail to reproduce the scaling violations diplayed 
by SIDIS data reported by the HI collaboration [12]. 

On year 2007, two parametrizations (HKNS [26] and DSS [27, 28], including also [29]), 
followed by AKK08 [30] one year later, have been released which include also data from 
SIDIS and hadronic collisions, and show an error analysis in the fit. Their main features 
are listed in Tab. I. 

The scaling violations are described by solving evolution equations at NLO with 
different techniques and constraining them to reproduce the longitudinal momentum sum 
rule. SU(2) isospin symmetry is assumed for the unfavoured (sea) channel at the starting 
scale; AKK08 and HKNS assume it also for the favoured one. AKK08 further includes 
the resummation of leading (LL) and next-to-leading (NLL) logarithms for z — >■ 1 both 
in the evolution equations and in the coefficient functions of the factorization formula, 
somewhat confusing the comparison with the other fixed-order extractions. 

In fact, AKK08 and DSS produce very similar results but at large z [27], where the 
effect of large logarithms is dominant; HKNS gives a doubtful softer gluon D'^^ because 
it lacks the constraint from RHIG pp data [26] . Remarkably, all sets fail to reproduce the 
STAR data for A, A production in pp collisions by almost one order of magnitude [23] . 

2'2. Tranverse-momentum, dependence. - The completely unknown dependence upon 
the transverse momentum of partons is usually parametrized in terms of a Gaussian 
ansatz. In SIDIS, the Gaussian width is fixed, for example, by reproducing the data for 
the average transverse momentum squared (P^j^) of final hadron h with respect to the 
virtual photon direction in the lab [31], because 



(1) 



(PL)=^'(pi) + (K|), KT = -^kT, 
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Ref : <k|> = 0.28, <pj> = 0.25 

h*: <k|> = 0.15 ± 0.02, <pj> = 0.45 ± 0.01 

h": <kj> = 0.06 ± 0.02, <p^> = 0.48 ± 0.01 



« 0.6 - 




COMPASS 2004 LiD (part) Preliminary 



0.4 - 




0.6 



0.7 

Z2 



Fig. 1. - The COMPASS 2004 data for {PIj_) as a function of z'^ [33]. Lower (higher) curves 
refer to fits for negative (positive) charged final hadrons. Straight (green) fine obtained with 
constant Gaussian widths. 



where p_L refers to the initial parton considered again in the lab frame, while is the 
hadron transverse momentum with respect to the direction of the fragmenting parton 
(approximately, the jet axis) and, viceversa, refers to the fragmenting parton with 
respect to the final hadron. 

A new combined analysis of recent SIDIS data (including azimuthally asymmetric 
cos (f> and cos 2(f) modulations in the cross section) and Drell-Yan data has lead to a 
new parametrization (see Ref. [32] and references therein). The best fit gives (p^) = 
0.38±0.06, (K^) = 0.16±0.01 GeV^, with a linear dependence in the cm energy s which 
broadens the distributions for increasing energy. 

However, there are several indications in SIDIS measurements that the Gaussian 
widths should depend at least on z and on the flavor content of the final hadrons. For 
example, in Fig. 1 the COMPASS 2004 data for (P^_l) [33] show a clear dependence on 

and on the hadron charge, the straight line being obtained with constant (p^), (K|,), 
and fitting only few data points. 

Similarly, in Fig. 2 the asymmetry between deuteron and proton targets of multiplic- 
ities for pions and kaons are displayed as functions of z for recent HERMES data [34]. 
Nonvanishing (sometimes, large) asymmetries indicate the marked sensitivity of final- 
state distributions to the different flavor content of different targets. 

Finally, in Ref. [35] a first attempt in describing the evolution of transverse-momentum 
dependent (TMD) nonperturbative functions, either initial distributions or final fragmen- 
tations, was put forward in the context of a proper factorization theorem. It turns out 
that even at LO, for the test case -Df „ of interest here, the Kt distribution strongly de- 
pends on the hard scale even at very low values of |Kt|, getting broader and broader 
with increasing Q^. 

2'3. Polarized fragmentation. - There is only one polarized fragmentation function 
that has been extracted from experimental data so far: it is the socallcd Collins function 
[2]. It is related to the probability density of having a distorted P;ij_ distribution 
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Fig. 2. - The HERMES data for the muhiphcity asymmetry between deuteron and proton 
targets as a function of z [34]. Higher (lower) panels for pions (kaons), left (right) columns for 
positive (negative) charges. 



of the final hadron h depending on the direction of the transverse polarization Sg of the 
fragmenting quark via the spin-orbit effect Sg • k x P/j, with k pointing in the direction 
of the jet axis. 

It can appear in e,^ e" — >■ q' — > h'^h^X events or, more interestingly, in SIDIS on 
transversely polarized targets. In fact, a specific azimuthally asymmetric modulation of 
the leading- twist SIDIS cross section contains the convolution h\®Hi J^ on the transverse 
momenta of the initial and final quarks, where hi is the socalled transversity parton dis- 
tribution, a poorly known cornerstone in the construction of the (spin) partonic structure 
of the nucleon (for a review, see for example Ref. [36]). 

Both azimuthal asymmetries in e^e~ and in SIDIS have been measured by the 
BELLE [37] and HERMES and COMPASS [38, 39] collaborations, respectively, and 
also later with increased statistics [40, 41, 42]. A simultaneous fit of the three data 
sets made it possible for the first time to extract a parametrization for the transversity 
hi [43, 44]. The results are presently limited to the valence u,d, quarks because of the 
limited kinematical range covered by the experiments. As for the fragmentation, we 
speak of favoured (u — > tt^, d — > tt^) and unfavoured {u — > tt^, d — > 7r+) channels with 
the surprising and interesting findings that Hi^ « —H^^^ [38]. 

As prevously said, the Collins function is paired to the transversity in a convolution 
on quark transverse momenta. Hence, the knowledge of the whole H^{z,K.t) function 
is crucial to unravel the convolution. However, from the expression of the azimuthal 
asymmetry in e+e^ only a K-r-integrated moment of can be isolated, leaving its 
Kt dependence unconstrained. The latter has been parametrized similarly to Z?i, i.e. 
with a Gaussian ansatz independent of kinematics and of the involved fiavor. Therefore, 
the extraction of hi in Ref. [43, 44] is affected by a model dependence. Moreover, the 
azimuthal asymmetries in e~^e~ and in SIDIS were measured at two very different scales, 

= 100 and 2.5 GeV^ respectively, but evolution effects in the Kt dependence were 
neglected. Partial results on the evolution properties of the Collins function have been 
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Fig. 3. - The spectator approximation for a parton with momentum k fragmenting into a detected 
hadron with momentum Ph- 



recently published [45], but it is fair to say that a full treatment of TMD evolution in 
the Collins effect is still missing. 

3. Models 

Since the extraction of fragmentation functions from experimental data is affected 
by large uncertainties, as we have seen about the Collins function and, more generally, 
about the dependence acquired by hadrons during the fragmentation, it is desirable 
that this phenomenology is supported by model speculations. In the following, we sketch 
three main classes of models that appeared in the recent literature. 

3'1. Spectator approximation. ~ The spectator approximation amounts to describe the 
fragmentation as the decay of a parton with momentum k into the observed hadron h 
with momentum P/j leaving a residual system in an on-shell state with momentum k — 
(see the diagram in Fig. 3). The latter condition grants that most of the calculations 
can be performed analytically, including the expression for the off-shellness fc^(z) of the 
fragmenting parton. The drawback is that only the favoured channel can be taken into 
account. 

For the typical u — 7r+ channel, two main choices have been adopted in the literature 
for the quark-pion-spectator vertex: the pseudoscalar coupHng g,rg75 [46, 47, 48, 49, 50] 
and the pseudovector coupling (?7rij757p-PC [51, 52, 48]. In all cases the coupling was 
assumed to be point-like except in Refs. [50, 49], where a gaussian form factor was used 
with a z-dependcnt cut-off. 

Complicated objects like the Collins function appear if there are nonvanishing in- 
terference diagrams involving different channels. In the spectator approximation, these 
final-state interactions can be achieved by adding to the left or right side of the diagram 
in Fig. 3 insertions involving pions and/or gluons. As an example, in Fig. 4 the K^- inte- 
grated i-moment H^^ ^^^'^^ (normalized to Df\^) from Ref. [49] is plotted as a function 
of z for three different hard scales and compared with the parametrization of Ref. [43] , 
whose statistical error is represented by the uncertainty band. The spectator results were 
obtained using a pseudoscalar qir coupling and gluon insertions. The model parameters 
were fixed by reproducing the unpolarized Di at the lowest available — 0.4 GeV^, 
as it was extracted from e+e^ data in Ref. [53]. Since the parametrization of Hi was 
performed using SIDIS data for the Collins effect at = 2.5 GeV^, the band in Fig. 4 
should be compared with the dashed (green) line, showing a substantial agreement with 
the spectator model. 
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Fig. 4. - The spectator model result for , 
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{z)/Dl^{z) from Ref. [49] (see text). Solid 



(black), dashed (green), and dot-dashed (red) lines for = 0.4,2.4,110 GeV^ 
Uncertainty band for the phenomenological extraction of Ref. [43]. 



respectively. 



3'2. Namhu-Jona Lasinio jet model. - In the Nambu-Jona Lasinio (NJL) jet model [54] , 
the fragmentation is represented as a recursive process depicted in Fig. 5. The D™^{z) 
can be obtained by solving a set of coupled integral equations based on the following 
product ansatz: 



(2) (z)dz = d^{z)dz + Y.[d'^^ D^] {z) , 

Q 

where the elementary fragmentation function d™ of a quark q in the meson m describes 
each splitting step in Fig. 5, and the sum runs upon all allowed intermediate states Q in 
the cascade. In each splitting, the depends on the quark-meson coupling gqmQ, which 
is determined from the residue (at the pole of the meson mass) in the quark-antiquark 
T matrix [55]. 

However, the above framework is justified only in the Bjorken limit where the quark 
initiating the cascade has an infinite momentum and produces an infinite number of 
hadrons; only in this limit the momentum sum rule is satisfied. Moreover, the general- 
ization to fragmentation into baryons is not trivial. Finally, solving the coupled integral 
equations is sometimes a heavy computational task. For all these reasons, the quark- 
cascade description of the fragmentation has been approached using the Monte Carlo 
technique [56]. In this context, the fragmentation function D'^{z)Az for a hadron h with 



^ ^ ^ 
• * • • • »— 

q a Q' 0" 



Fig. 5. - Recursive quark-meson splitting as the elementary mechanism of fragmentation in the 
Nambu-Jona Lasinio jet model. 
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Fig. 6. - The zD*l^^{z) for h — tt"*" (left panel) and h = tt~ (right panel). Dashed and dot-dashed 
lines for the parametrizations of HKNS [26] and DSS [27], the uncertainty bands referring to 
HKNS. Dotted line for the result of the NJL jet model, solid (red) line includes in the cascade 
also intermediate vector mesons [56]. 



momentum fraction in the range [z, z+Az], is deduced by calculating the average number 
of hadrons of type h produced in the cascade depicted in Fig. 5 for a predefined number 
of steps iVstop- Each step of the cascade is randomly sampled using the calculated 
in the NJL jet model, and the entire cascade is simulated Nsim times with Ngim large 
enough to stabilize the average. When the hadron /i is a baryon, the internediate states 
Q are described in the framework of the spectator scalar diquark model [56] . 

In Fig. 6, the zDi^^{z) is shown for the favoured {h — 7r+, left panel) and unfavoured 
channels {h — 7r~, right panel). The dashed and dot-dashed lines represent the phe- 
nomenological parametrizations of HKNS [26] and DSS [27] (see Tab. I), the uncer- 
tainty bands referring to HKNS. The result of the NJL jet model is represented by the 
dotted line, while the solid (red) one includes in the cascade also intermediate vector 
mesons, further decaying into the observed tt*. For sake of consistency with the em- 
pirical parametrizations, the NJL jet model results are evolved at NLO from Qq — 0.2 
to = 4 GeV^. The Lepage-Brodsky scheme adopted in regularizing the calculation 
of loop diagrams reduces the available range in z depending on the hadron type h. For 
z Zynin{h) the Monte Carlo result for zDi{z) tends to a constant which becomes larger 
for increasing Q^, from which the apparently divergent behaviour for small z. We refer 
the interested reader to Ref. [56] for the results about other channels. 

3'3. Recursive model. - Present Monte Carlo event generators of quark and gluon 
jets do not include spin in the elementary degrees of freedom. Therefore, in order to 
have a guide in studying azimuthal asymmetries like the Collins eifect, a new quantum 
approach to polarized quark fragmentation was suggested in Ref. [57]. As an example, the 
production of N pseudo-scalar mesons /ii,/i2, . . . ,/iAr, in e'^e~ annihilation is depicted 
in Fig. 7 via the elementary annihilation qoq^i, where in each step i of the chain, from 
right to left, we have qi-i — hi + qi for i = I, .., N, and qN = —q-i- 

Considering the upper part of Fig. 7 as a Feynman amplitude with spinors, vertices 
and propagators, it is possible to study its dependence on spin by introducing some ad hoc 
simplifications. In Ref. [57], this amplitude was estimated in a socalled multiperipheral 
model by using Pauli spinors and matrices in the vertices, and by approximating each 
intermediate fermion propagator at step i with an expression similar to the meson-nucleon 



UNPOLARIZED AND POLARIZED FRAGMENTATION FUNCTIONS 



9 




Fig. 7. - The process e^e — > go9-i — > /t-i + '^«2 + • ■ ■ + ftjv as a recursive q — > hq' splitting. 



scattering amplitude, 

(3) Ai « exp[-6h2j,/2] [^(h^j,) + id • z x h.^] , 

i.e. with a non-spin-flip complex function fj, and a spin-flip part, b being some free 
parameter. These prescriptions can be shown to respect invariance under all "good" 
transformations like rotations, boosts, and parity, all considered with respect to the jet 
axis z. 

If Im(/i) 7^ 0, this imaginary part can be shown to act as a source of transverse 
polarization at step i even if the quark was unpolarized or longitudinally polarized at 
step i — 1 [57] . This means also that during the cascade the helicity of a quark can be 
partly converted to its transversity or viceversa. As a consequence, if Im(/i) ^ one can 
have for = 1 a Collins effect Si • z x hiT, and for = 2 an iterated Collins effect 
with alternate sign, which could explain the experimental flndings H^^"^^ ss —H^^^^ 
described in Sec. 2'3 [38]. This result confirms the outcome of the Lund ■^Pq string 
mechanism [58] . But in addition it contains the three-particle correlation z • h2T x hi^ 
named jet handedness [59], which is interpreted as a two-step mechanism: at z = 1, a 
transverse polarization Si^ || hiT is generated from the helicity Sqz of previous step; at 
i = 2, a Collins effect takes place as z-h2T x Sit, which coincides with the jet handedness. 

Further work is needed to promote the multiperipheral model of Ref. [57] to a real- 
istic Monte Carlo event generator. For example, one should include antiquarks in the 
fragmentation cascade, or explore the interference of the amplitude in Fig. 7 with dia- 
grams showing differently ordered A^ hadrons. Preliminary experimental results already 
appeared for SIDIS production by the HERMES collaboration (an almost vanishing 
Collins effect [60] and a large cos 2;/) asymmetry in the unpolarized cross section [61]) 
that cannot be easily accommodated in the multiperipheral model in its present version. 

4. Di-hadron Fragmentation Functions 

As already sketched in Sec. 2'3, the extraction of the transversity parton distribution 
via the Collins effect suffers from several uncertainties and model dependencies, mostly 
related to the need of dealing with TMD objects. A complementary approach is provided 
by the semi- inclusive process ep^ e'{hih2)X where two unpolarized hadrons with 
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momenta Pi and P2 emerge from the fragmentation of the same quark. The kinematics 
is similar to the single-hadron SIDIS except for the final state, where the hadron pair 
carries a fractional energy z — Z1+Z2 with a total momentum P/j = P1+P2 and a relative 
momentmn R ~ {Pi— P2) /2 (see Fig. 8, left panel). The underlying mechanism is related 
to • X R: the transverse polarization of the fragmenting quark q is transferred to 
the relative orbital angular momentum of the hadron pair. Contrary to the Collins 
effect, this mechanism survives after integrating away the transverse momentum of each 
particle and can be analyzed in the coUinear factorization scheme [3]. The probabilistic 
weight for this to happen is represented by the polarized dihadron fragmentation function 
Hf^^^^{z,M'^,Q'^), where Pf^ = ^ is the pair invariant mass and represents a 
new soft scale in the process. It is easy to show that |Rp is a linear function of [63]. 

Dihadron Fragmentation Functions (DiFF) were introduced in Ref. [64] and studied 
for the polarized case in Refs. [65, 66, 67]. The decomposition of the SIDIS cross section 
in terms of parton distributions and DiFF was carried out at leading twist in Ref. [68] 
and to sub-leading twist in Ref. [69]. Since Rt = Rsinfl, where in cm frame of the 
hadron pair 6 is the angle between Pi and the direction of Ph in the lab frame [63], the 
leading-twist cross section shows an azimuthally asymmetric modulation proportional to 
sin{(j)fi + sin 0, where 4>fj, 4>Si are defined in Fig. 8. The proportionality coefficient 
contains the product hf H^^^^'^^ [66, 3, 63, 68, 70]: the advantage of working in coUinear 
factorization scheme reflects in a very simple relation with no convolution on transverse 
momenta, as in the case of the Collins effect. 

The H^^'^^'^^ is sensitive to the interference between the fragmentation amplitudes into 
hadron pairs in relative s wave and in relative p wave [63] . Intuitively, if the fragmenting 
quark is moving along z and is polarized along y, a positive H^^'^''^^ means that hi is 
preferentially emitted along — x and /12 along x. The corresponding unpolarized partner 
D'I^^^ is averaged over quark polarization and hadron pair orientation. Similarly to the 
single-hadron SIDIS, the unknown DiFF must be independently determined from e~''e~ 
annihilation producing, in this case, two hadron pairs (see right panel of Fig. 8, adapted 
from Ref. [62] ) . The relevant signal is similar to that of the Collins function, except that 
each transverse polarization of the qq pair is now correlated to the azimuthal orienta- 
tions (t>ii, 4>]i, of the planes formed by the momenta of the corresponding hadron pairs. 
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Fig. 9. - The combination of valence u, d, flavors for transversity from Eq. 4. Black circles 
for the SIDIS data from HERMES [72], red squares from COMPASS [77]. The error bars are 
obtained by propagating the statistical errors of each term in the equation. The uncertainty 
band represents the same observable as deduced from the parametrization of Ref. [44] . 



suggesting that is related to the above mentioned concept of jet handedness [59, 71]. 
In the leading-twist cross section, this correlation shows up as a modulation proportional 
to cos((/)fl, + (jjii) [71]. 

The sine ^^^^ asymmetry in SIDIS production of a tt+tt" pair off a trans- 

versely polarized proton target was measured for the first time by the HERMES collabo- 
ration [72] , ruling out the model of Ref. [66] and showing compatibility with predictions 
based on the spectator approximation [73, 74] (see also the later Ref. [75]). Prelim- 
inary results are available also from the COMPASS collaboration using deuteron [76] 
and proton [77] targets. The ^'=°s(0r+</>j5) asymmetry in e+e~ annihilation was recently 
measured by the BELLE collaboration [62]. By combining these data, the transversity 
hi was extracted for the first time in a coUinear factorization scheme, properly including 
LO evolution effects in the behaviour of DiFF at different scales [78] . 

The result is shown in Fig. 9. By assuming certain symmetries of DiFF for tt+tt^ pairs 
under isospin and charge conjugation transformations [73, 78], the SIDIS spin asymmetry 
^sin(0R-(-0s) sine ^^^^ written, after integrating upon z,M^, as 



K-{x,Q')--hi^{x,Q^) 



(4) 



A 



UT 



E 



where x is the fractional momentum carried by quarks, Cy is the depolarization factor, 



nu{Q^)= / dzdMlDlJ~{z,MlQ'') 
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(5) n^(Q^) = I dzdM'J^^Ht/--{z,MlQ') . 

The densitiy (n„) of tt+tt" pairs produced by a (un)polarized quark u can be calcu- 
lated by fitting the (z, M^) dependence of DiFF in the e+e~ asymmetry as 
described in Ref. [78]. The integrals in Eq. 5 are performed over the range 0.5 < M/j < 1 
GeV, 0.2 < ^ < 0.7, considered by HERMES [72]. Evolution effects at LO [79] produce 
a reduction of by 92%± 8% when moving from the BELLE scale = 100 GeV^ 

down to the HERMES (Q^) = 2.5 GeV^. Using the MSTWLO08 set of unpolarized par- 
ton distributions /i [80], the data points of Fig. 9 are obtained by inserting experimental 
data in ^^^("^H+flis) sme q-j- ^ Black circles correspond to the HERMES measurement 
of Ref. [72], red squares to the COMPASS one of Ref. [77]. The error bars are obtained 
by propagating the statistical errors of each term in Eq. 4, the major role being played by 
the experimental statistical errors on ^^^('?^K+'?*s)sm6i rpj^^ uncertainty band represents 
the same observable as deduced from the parametrization of Ref. [44] . From Fig. 9 we 
deduce that there is a substantial agreement between the two extractions of hi obtained 
from two independent methods. However, the statistical meaning of the displayed error 
bars is totally different from the uncertainty band of Ref. [44]: in order to perform a 
meaningful comparison a more detailed analysis is needed [78]. 

5. — Outlooks 

There are several interesting ongoing developments in each of the fields touched in 
previous Sections. 

As for single-hadron fragmentations, we are rapidly moving towards a full NNLO anal- 
ysis of evolution effects when connecting data for Z?^ j(z, Q^) at different scales [81, 82]. 
For h = K, nonsinglet fragmentation functions for can be directly extracted from 
data in a model independent way [83] , the present limitation being due to weak constrains 
coming from a not enough large data sc^t. The ultimate goal would be to construct 
a nonperturbative error for D'l^ by comparing in a meaningful way the most reliable 
parametrizations (like HKNS, DSS and AKK08, described in Tab. I), similarly to what 
is done for parton distributions. To do so, a common interface similar to LHAPDF is 
needed, while at present only the web site http : //www . pv . inf n . it/~radici/FFdatabase/ 
is available. The dependence of D^^ upon the transverse momentum of h is probably 
the most developing field, since very few experimental inputs are available. New results 
have been recently obtained which open the door to a correct treatment of evolution 
effects for TMD nonperturbative soft functions in the context of a suitable factorization 
theorem [35]. 

As for two-hadron fragmentations, first data from the BELLE collaboration were 
released very recently and much more work has to done in order to unravel dihadron 
fragmentation functions and to confirm the first results about the extraction of transver- 
sity. In particular, a full flavor separation is in order for the analysis of e+e" — >■ 
(7r+7r~)(7r+7r~)X data, enlarging the explored range in invariant mass of the pion pair 
to study more contributing resonances. Moreover, new data from the COMPASS col- 
laboration have been released [84], which demand for a more refined analysis. Finally, 
some data were released from the PHENIX collaboration about the pp^ (7r+7r^)X 
process [85] (and more data will probably become available from the STAR collabora- 
tion in the near future), which should help in separating the antiquark components of 
dihadron fragmentation functions and, consequently, of transversity [86] . 
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